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\D '. Abstract 

The effects of u> meson self-coupling (OMSC) on the properties of finite nuclei and neutron 
(3J[)' stars are investigated within the framework of effective field theory based relativistic mean-field 

■ (ERMF) model which includes the contributions from all possible mixed interactions between the 

scalar- isoscalar (er), vector-isoscalar (u) and vector- isovector (p) mesons upto the quartic order. 



X 



For a realistic investigation, several parameter sets corresponding to different values of OMSC 
are generated by adjusting the remaining parameters of the ERMF model to fit the properties of 



> 

the finite nuclei. Though, all these parameter sets give equally good fit to the properties of the 
finite nuclei, only moderate values of OMSC are favored from the "naturalness" point of view. The 
equation of state for the symmetric nuclear and pure neutron matters resulting from the parameter 
sets with the moderate values of OMSC are in close agreement with the ones obtained within the 
Dirac-Brueckner-Hartree-Fock approximation. For such parameter sets the limiting mass for the 
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neutron stars composed of /3-stable matter is ~ 1.9M@. It is found that the direct Urea process 
can occur in the neutron stars with "canonical" mass of 1AM @ only for the moderate and higher 
values of OMSC. Some other interesting properties for the neutron stars are also discussed. 

PACS numbers: 21.10.-k,21.65+f,24.30.Cz,21.60jz,26.60.+c 
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I. INTRODUCTION 



The concepts of effective field theory (EFT) have provided a modern perspective to the 
relativistic mean-field models (jj, |2j, |3|. The EFT based relativistic mean-field (ERMF) 



models are obtained by expanding the energy density functional in powers of the fields for 
scalar-isoscalar (a), vector-isoscalar (u) and vector- isovector (p) mesons and their derivatives 
upto a given order v . Thus, a ERMF model includes the contributions from all possible self 
and mixed interaction terms for a, uo and p mesons in addition to the cubic and quartic self 
interaction terms for a meson aspresent in the conventional quantum hydrodynamic based 
relativistic mean field models 4], [J. The parameters (or the expansion coefficients) appearing 
in the energy density functional are so adjusted that the ERMF results for a set of nuclear 
observables agree well with the corresponding experimental data. The EFT demands that 
the adjusted parameters must exhibit "naturalness", i.e., the values of all the parameters 
should be roughly of same size when expressed in appropriate dimensionless ratios 6]. The 
lack of "naturalness" implies that the omitted terms are important. Sometimes not all the 
terms upto a given order v are considered which might lead to unnaturalness [l]. It is found 
in Ref. j^] that the ERMF models containing terms upto order v = 4 can be satisfactorily 
applied to study the properties of finite nuclei. The inclusion of next higher order terms 
improves the fit to the finite nuclear properties only marginally. In Ref. (f| it has been shown 
that even high density behaviour of the equation of state (EOS) of pure neutron matter and 
/3-stable matter are predominantly controlled by the quartic order uo meson self-coupling 
(OMSC). The effects of inclusion of higher order terms on the high density behaviour of 
EOS are found to be only modest to negligible. 

There have been several parameterizations of ERMF models containing most of 

the terms upto the quartic order [y = 4). These parameter sets are obtained by a fit to 
the set of experimental data for the properties of a few stable closed shell nuclei. The 
fitted parameters exhibit "naturalness". But, some of the nuclear matter properties they 
yield requires attention; the linear density dependence of the symmetry energy coefficient is 
strong and the nuclear matter incompressibility coefficient K is either little too low or quite 
high. The linear density dependence of the symmetry energy coefficient and the nuclear 
matter incompressibility coefficient should be adjusted to yield reasonably the neutron-skin 
thickness and the centroid energy of the isoscalar giant monopole resonance, respectively. 
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The FSUGold parameter set is obtained recently 0, |^| in such a way that they give realistic 
values for all the properties normally associated with the nuclear matter at the saturation 
density. However, some of the parameters of FSUGold set show deviation from "natural- 
ness" . This may be pointing to the fact that the contributions due to the mixed interaction 
terms involving o — uo and o — p mesons omitted from the FSUGold parametrization are 
important. For FSUGold parameters, the limiting mass of the neutron star is only 1.7 Mq 
which might reduce further with the inclusion of hyperons degree of freedom |ld. lll|. On 
the other hand, an early study using ERMF model shows Ref . |6[ that by changing the value 
of OMSC within a reasonable range, the limiting mass of neutron stars can be varied from 
1.8M - 2.8M for a fixed set of nuclear matter properties. Of course, for a realistic variation 
in the value of limiting mass, the nuclear matter properties must be determined from the 
experimental data on the finite nuclei. It is also to be noted that in Ref. [f| the effects of 
the mixed interactions were not considered and no constrain was imposed on the density 
dependence of symmetry energy coefficient. 

In the present work we have performed a realistic investigation involving the effects of 
OMSC on the properties of finite nuclei and neutron stars using a ERMF model. We 
consider a ERMF model which includes the contributions from all possible mixed interactions 
between the ex, u and p mesons upto the quartic order. Like, in the case of FSUGold 
parametrization, we determine all the properties normally associated with nuclear matter 
from the experimental data on the finite nuclei. We generate several parameter sets for fixed 
values of OMSC by adjusting the remaining parameters of the ERMF model to fit exactly 
the same set of experimental data for the total binding energies and charge rms radii for some 
closed shell nuclei. The binding energy data is considered for nuclei ranging from normal to 
the exotic ones. We also include in our fit the value of the neutron-skin thickness for 208 Pb 
nucleus to constrain the linear density dependence of the symmetry energy coefficient. We 
restrict the value of the nuclear matter incompressibility coefficient to be within 220 — 240 
MeV as required by the experimental data on the centroid energies for the isoscalar giant 
monopole resonance. The best fit parameters are searched using the simulated annealing 
method (SAM) which we have applied recently ^ (3 to determine the parameters of the 
standard and generalized Skyrme type effective forces. 

In Sec. |H]we discuss ERMF model in brief. In Sec. II I II we discuss the simulated annealing 
method used in the present work to search for the best fit parameters. Several parameter 
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sets for the different values of the OMSC generated in the present work are given in Sec. 
IIVI In Sees. IS - 1 VI II we discuss our results for the symmetric nuclear matter, pure neutron 
matter, finite nuclei and neutron stars obtained using our newly generated parameter sets 
and compare them with the ones for the FSUGold parametrization. As a customary, we 
also compare our results with those obtained using most popularly used NL3 parameter set 
[l^ . Finally, in Sec. IVIIII we summarize our main conclusions. 

II. ENERGY DENSITY FUNCTIONAL 

The derivations of the EFT motivated effective Lagrangian and corresponding energy 
density functionals are weU do _ d ln Refs. fl QU m tMs section we shaU out^e 
very briefly the final expressions for the energy density functional and field equations within 
the mean-field approximations as used in our numerical calculations. The energy density 
functional containing all the self and mixed interaction terms for the a, u and p mesons 
apto the ^ico^ea^Ea, 

£{r) = {-«* ■ V + (3 [M - $(r)] + W(r) + \r z R{r) + I±^A(r) J M (r) 



+ i^$ 2 (r) + -$ 3 (r) + — $ 4 (r) - ^W\r) - —R A (r) + J-[V$(r)f 
2 g 2 a v ; 6 v ; 24 K J 24 v ; 24 2g%. v n 
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~ ^[V^Mf - -^-W 2 {r) -a^{r)W 2 {r) - -a'^ 2 (r)W 2 (r) - --£R 2 (r 

- -L[Vi?(r)] 2 + a 2 $(r)i? 2 (r) - \a' 2 <$> 2 {r)R 2 {r) - \a^W 2 {r)R 2 {r) 
^9p 2 2 



2^[VA(r)] 2 



where, the index p runs over all occupied states of the positive energy spectrum. The 
variables $, VFand R represent the a, u and p meson fields, respectively, and A is the 
photon field. We must emphasize at this point that in the present work we consider all the 
mixed interaction terms. 

The mean-field equations for the nucleons, mesons and photons are as follows, 



-ia • V + (3 [M - $(r)] + W(r) + ^r 3 R(r) + ^pA(r) } o„( r) = r /( o /( (r). (2) 
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K A 

-A$(r) +m^(r) = g 2 a p s (r) - -$ 2 (r) - -$ 3 (r) + a^ir) 

I o 

+ a\<S>(r)W 2 (r) + a 2 R 2 (r) + a 2 <S>(r)R 2 (r] 



r 

-AW(r) + m 2 W(r) = g 2 p(r) - j:W 3 {r) - 2aMr)W(r) 
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(3) 



- a' 1 $ 2 (r)iy(r) - a' 3 W(r)R 2 (r 



-A/?(r) m 2 p R(r) = ^g 2 p p 3 (r) - |i? 3 (r) 27T 2 <I>( /•)/?(/•) 
-a' 2 $ 2 (r)i?(r) - a' 3 R(r)W 2 (r) , 



(4) 



(5) 



-AA(r) = e 2 p p (r). (6) 

The different kinds of baryon densities p s (scalar), p (total), p 3 (isovector) and p p (proton) 
appearing in Eqs. © - © are given by, 

P.(r) = X>£(r)/ty M (r), (7) 

p(0 = X>iM^M, (8) 

P 3 (r) = X>^)T 3 <fV(r), ( 9 ) 
n 

The EFT imposes the condition of "naturalness" on the parameters or the expansion coef- 
ficients appearing in Eq. ((TJ. The "naturalness" implies that the coefficients of the various 
terms in the energy density functional, when expressed in appropriate dimensionless ratios, 
should all be of the same size. The dimensionless ratios are obtained by dividing the Eq. 

by M 4 and expressing each term in powers of $/M, W/M and 2R/M 
that the dimensionless ratios , namely, 

1 1 1 k A 



0. 



This means 



2c 2 M 2 ' 2c 2 M 2 ' %c 2 n M 2 ' 6M' 24' 

(11) 

£ £ Q!i «2 «2 <^3 

24' 384' M' T' 4M' T' ~8~ 
should be roughly of same size. In the above equation c 2 = g 2 /m 2 , where i represents a, u 

and p mesons. 
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III. SAM ALGORITHM FOR X 2 MINIMIZATION 



The SAM is analogous to an annealing process in which a metal, initially at high tem- 
perature and disordered, slowly cools so that the system at any time is in a thermodynamic 
equilibrium. As cooling proceeds, the system becomes more ordered and approaches a frozen 
ground state at zero temperature. It is an elegant technique for optimization problems of 
large scale, in particular, where a desired global extremum is hidden among many local 
extrema. The SAM has been found to be an extremely useful tool for a wide variety of 
minimization problems of large non- linear systems in many different areas of science (e.g., 
see Refs. jisl 13, Q ) . Recently 18, [3], the SAM was used to generate some initial trial 
parameter sets for the point coupling variant of the relativistic mean-field model. 

In the present work we use SAM to minimize the \ 2 function defined as, 

1 Nd / M exp M th \ 2 



x 2 



where, Nd is the number of experimental data points and N p the number of fitted param- 
eters. The Oi stands for theoretical error and M^ xp and Mf are the experimental and the 
corresponding theoretical values, respectively, for a given observable. Since, the Mf 1 in 
Eq.flT!Z|) is calculated using the ERMF model, the values of \ 2 depends on the values of the 
parameters appearing in Eq. In Ref.[l2| we have described in stepwise manner the 

implementation of the SAM algorithm to determine the best fit parameters for the standard 
Skyrme type effective forces. To apply SAM one needs to specify the range for each of the 



parameters needs to be fitted and a set of guess parameters. In Ref. [12| the range of the 
Skyrme parameters were specified in terms of the properties of the nuclear matter at the 
saturation density. This was possible, because, the number of the Skyrme parameters to be 
fitted were directly related to the equal number of nuclear matter properties which led to a 
drastic reduction in the parameter space so that the best fit parameters could be searched ef- 
ficiently. For the case of generalized Skyrme type effective forces we used the SAM algorithm 
bv specifying the parameter space directly in terms of the range of each of the parameters 
13j . Since, the number of parameters for the generalized Skyrme type effective forces are 
larger than the number of the quantities normally associated with the nuclear matter. For 
the ERMF model we have used rather a hybrid approach. The parameters g a , g w , k and A 
of Eq.tJTJ) for a given values of a\, a{ and ( are expressed in terms of the binding energy per 



nucleoli e, incompressibility coefficient K, effective mass M* and saturation density po for 
the nuclear matter as follows jl , 20 ] , 



a 



2^ 



1 \§ 2 U" - 3W + 6U 



(13) 



C 2 



w 



(14) 



-$ 2 [7" + 5QU' - 8U 

I$3 



(15) 



where, 



_ \<$> 2 U" - 2&U' + 3U 

24 ^ 



$ = M — M* 



(16) 



(17) 



W = e + M- a k 2 +M* 2 



(18) 



(19) 



with 



C/" = p' s -( 2 1 , + (5' 1 H' 2 + 



(f 



2fc F _B* ~ u 9p 



/ = _2W («1 + «!$) 



(20) 



(21) 



1 



(22) 



-± E + {2a l + a' 1 <$> 2 ) + \(W 2 
In Eqs. (fT3|) - (J22J), all the quantities are obtained at the nuclear matter saturation density 

P = Pa- 
in Table U] we specify the parameter space in terms of e, K, M* , p and some other 
parameters. It can be easily verified using second and third columns of Table H] that the 
range of each of the parameters a±, a[, a?2, «2 an d «3 becomes identical (0.0 — 2.0 x 10~ 3 ) 
when expressed in appropriate dimensionless ratios using Eq. ljlljl . The column labeled "d" 
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denotes the maximum change allowed in a single step for a randomly selected parameter. 

The guess values of the parameters which we have used for initiating the SAM are given in 

last column of Table UJ The values of guess parameters are so chosen that they lie in the 

middle of the specified range. It should be noted that the parameters ( and £ denoting the 

self couplings for u and p mesons do not appear in Table E as they are kept fixed during x 2 

minimization procedure . The best fit parameters for which % 2 is the minimum are searched 

within the specified parameter space with the help of annealing schedule. The annealing 

schedule determines the value of the control parameter T(k) at k— th step starting from its 

flfl 

initial value T(0). As in Refs. |12j, 112! we use Cauchy annealing schedule given by, 

T( k) = m (23) 

where, starting value of k is equal to zero and is increased in steps of unity after each 100N p 
reconfigurations or 10N p successful reconfigurations which ever occurs first. We keep on 
reducing the value of T using Eq. (j2Sj) until the efforts to reduce further the value of x 2 
becomes sufficiently discouraging. We have performed test calculations for two different 
values of T(0) taken to be 2.5 and 5.0. For both the test calculations the final value of 
the x 2 function is almost the same. In what follows, we present our results obtained using 
T(0) = 5.0. 



IV. NEW PARAMETRIZATIONS OF THE ERMF MODEL 

We search for the best fit parameters of the ERMF models by following the x 2 mini- 
mization procedure as briefly outlined in Sec. 11111 Our data set contains the total binding 
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22, 



231 ]. We consider 



energies and charge rms radii for several nuclei taken from Refs 
total binding energies for 16 ' 24 0, 40 < 48 Ca, 56 < 78 Ni, 88 Sr, 90 Zr, ioo,ii6,i32g n and 2o 8pb nuclei; 
charge rms radii for 16 0, 40 ' 48 Ca, 56 Ni, 88 Sr, 90 Zr, 116 Sn and 208 Pb nuclei . The theoretical 
error needed to evaluate the x 2 function ( Eq. IT2|l are taken to be 1.0 MeV for the total 
binding energies and 0.02 fm for the charge rms radii. In addition, we also fit the value 
of neutron-skin thickness for 208 Pb nucleus to constrain the linear density dependence of 
symmetry energy coefficient. The very accurate data of neutron-skin thicknesses are still 
not available. It is shown j^J in a relativistic mean-field based random phase approximation 
that the neutron-skin thickness of 0.175 fm in 208 Pb nucleus and K « 240 MeV are required 
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to adequately reproduce the centroid energies of the isoscalar giant monopole and isovector 
giant dipole resonances. A Dirac-Brueckner-Hartree-Fock (DBHF) calculation [25| using a 
realistic potential predicts the value of neutron-skin thickness in 208 Pb nucleus to be 0.188 
fm. However, recently extracted value of neutron-skin thickness for 208 Pb nucleus from the 
isospin diffusion data lie within 0.16 — 0.24fm indicating large uncertainties |26j. In our fit 
we use the value of neutron-skin thickness to be 0.18 ±0.01 fm for the 208 Pb nucleus. To this 
end, we must point out that for our parametrizations, the center of mass correction to the 
total binding energy E cm is evaluated within the harmonic oscillator approximation which 
gives E cm = \f\w and we take Hlu = 45 A~ 1/3 - 25A~ 2/3 MeV. 

The parameters ( and £ corresponding to self-couplings for uj and p mesons can not be 
very well constrained by the properties of finite nuclei. These parameters mainly determine 
the high density behaviour of the EOS . The impact of the parameter £ is found to be 
appreciable for pure neutron matter only at very high densities [6] . The seven new parameter 
sets FSUGZ00, FSUGZ01,....,FSUGZ06 are generated for fixed values of OMSC parameter 
C = 0.00, 0.01, 0.06 and keeping £ = 0. All the FSUGZ parametrizations give equally 
good fit to the properties of the finite nuclei. In Table |H] we give the values only for 
FSUGZ00, FSUGZ03 and FSUGZ06 parameter sets. In Table ITTT1 we present the values of 
these parameters when expressed in dimensionless ratios using Eq. (II II) . In Tables ITT1 and Hill 
we also list the values of the parameters for NL3 and FSUGold jjl^sets. It is clear from 
the Table ITTT1 that for the case of FSUGold parameter set the a' 3 is quite large in comparison 
to the other parameters when expressed in appropriate dimensionless ratios. We notice from 
Table ITTT1 that the values of the parameters 7c, A, 5i and a[ are strongly correlated with the 
values of (. For smaller (, the parameter A becomes small and negative which not only gives 
rise to deviations from the "naturalness" behaviour but also gives instabilities in the EOS 
at high densities. Whereas, for larger (, the values of 7c, a% and a[ become smaller and once 
again show deviations from the "naturalness". Thus, it appears that a moderate value of 
C ~ 0.03 is favored from the "naturalness" point of view. 

To understand better, the overall "naturalness" behaviour of the various parametrizations 
of the ERMF model presented in Table ITTT1 we consider the ratio of the largest to the smallest 
parameters for a given set. One expects this ratio to be of the order of unity if the parameters 
strictly obey the "naturalness" . We find that the ratio of the largest to the smallest param- 
eters are 9.8, 8.0 and 220.1 for the FSUGZ00, FSUGZ03 and FSUGZ06 parametrizations, 
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respectively. Though, the FSUGZ03 parameter set is favoured from the "naturalness" point 
of view, it requires further improvement to obey the "naturalness" criteria in a more strin- 
gent manner. From the naturalness viewpoint, FUGZ03 parametrization can be improved 
by considering few additional terms in the energy density functional (Eq.([Tj)). Indeed, it is 
clear from the findings of Ref. [3j that the next higher order term containing gradients of 
the fields are necessary for the overall improvements of the "naturalness" behaviour of the 
parameters. 

Before embarking on the discussion of our results, we would like to focus on the strong 
correlations existing between the parameters of the ERMF model which can be seen from 
Tables HU and HH It is evident from Eqs. flUID - 422) that for given values of e, K, M* and 
Po, the coupling constants g a , g^, k and A are correlated with ax, ct[ and (. The value of 
£ also affects the symmetry energy coefficient and its density dependence due to the mixed 
term containing u and p meson fields. In other words, the coupling constants g p and a' 3 
are also correlated with (. So, it is very important to constrain the value of (. This can 
be achieved only if the behaviour of the EOS at higher densities or the maximum mass of 
the neutron stars are known. To constrain the values of a\ and a[ one needs to know the 
density dependence of the effective meson masses. The values of a^, o.' 2 and a' 3 can be fixed 
by the neutron-skin thicknesses of asymmetric nuclei. 



V. NUCLEAR AND NEUTRON MATTERS 



In this section we discuss our results for the symmetric nuclear matter (SNM) and pure 
neutron matter (PNM). In Table HVl we present our results for the various properties asso- 
ciated with the SNM at the saturation density. The quantity L in Table II VI determines the 
linear density dependence of the symmetry energy coefficient J and is given by, 

dJ 



L = 3p 

dp 



(24) 



P=PO 

All the nuclear matter quantities given in the table for the FSUGZ and FSUGold parameters 
are more or less close to each other. The values of nuclear matter incompressibility coefficient, 
the symmetry energy coefficient and the linear density dependence of symmetry energy for 
the NL3 parameter set are much higher than that for our new parameter sets. 

We use our newly generated parameter sets to study the effects of OMSC parameter ( 



10 



on the EOS for the SNM and PNM. In Fig. ^we plot the difference between the energy per 
nucleon Ae = e(£ = 0) — e(£) as a function of ( for SNM(dotted lines) and PNM(solid lines) 
at fixed densities p = 0.5 fm -3 and 0.8 fm -3 . We calculate e(() at ( — 0.0, 0.01, 0.06 
using the parameter sets FSUGZ00, FSUGZ01,...., FSUGZ06, respectively. We see that Ae 
increases with increase in (. In other words, the e(£) decreases with increase in ( implying 
that EOS becomes softer as ( increases. At higher densities the values of e appears to be 
very much sensitive to the value of (. For example, the value of Ae at C — 0.06 is less 
than O.lMeV for p = 0.2 fm~ 3 and it becomes greater than 125 MeV for p = 0.8 fm -3 . It 
must be pointed once again, for different values of ( the remaining parameters of the ERMF 
models are so adjusted that they appropriately fit the properties of the finite nuclei. We find 
that FSUGZ03 parametrization which corresponds to a moderate value of ( yields an EOS 
which is neither too stiff nor too soft. In Figs. Ufa) and|^b) we plot the EOS for SNM and 
PNM obtained using FSUGZ03 parametrization which agree reasonably well with the EOS 
obtained within the DBHF framework j^. The NL3 and FSUGold parametrizations yield 
the EOS which show larger deviations, in particular, at higher densities as compare to the 
DBHF results. The NL3 parameters gives very stiff EOS while the FSUGold parameters 
gives rise to very soft EOS. 

VI. FINITE NUCLEI 

The fit to the properties of the finite nuclei for the parameter sets FSUGZ00 - FSUGZ06 
are very much similar. For the sake of clarity we present our results obtained for the finite 
nuclei only for the parameter set FSUGZ03. In Fig|3]we display the values of the relative 
error in the total binding energies SB = (B exv — B th )/B exp in percent obtained for the 
FSUGZ03 parametrization. For the comparison, we also display the similar results obtained 
for the NL3 and FSUGold parameter sets. It is clear that the binding energies obtained 
using FSUGZ03 parameter set agree better with the experimental data than those for the 
NL3 and FSUGold parameter sets. The rms errors in the total binding energies are 1.6, 2.5 
and 3.6 MeV calculated using FSUGZ03, NL3 and FSUGold parameter sets, respectively. 
For other FSUGZ parameter sets the rms errors in the total binding energies lie in between 
1.5 — 1.8 MeV. In Fig. H] we display our results for the relative error Sr c h in percent for 
the charge rms radii and compare them with the ones obtained for the NL3 and FSUGold 
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parametrizations. The values of charge rms radii calculated for various parametrizations 
differ only marginally. The rms errors in the charge rms radii are 0.03, 0.02 and 0.03 fm for 
the FSUGZ03, NL3 and FSUGold parameter sets, respectively. Our results for the neutron- 
skin thickness Ar = r n — r p are plotted in Fig. The values of the rms radii r n and r p 
are calculated using the point density distributions for neutrons and protons, respectively. 
The values of Ar shown by crosses are the predictions based upon the DBHF calculations 
(25! and isospin diffusion data [26j. We see that the values of Ar obtained using FSUGZ03 
and that for the FSUGold parameter sets are closer to the predicted values. The NL3 
parameters significantly overestimates the values of Ar. We would like to mention here that 
for the NL3 type of parametrization the values of the symmetry energy coefficient J and its 
linear density dependent L are determined by a single parameter g p which determines the 
coupling of p mesons with the nucleons. It may be therefore not possible with NL3 type of 
parametrization to fit simultaneously the values of the binding energies and the neutron-skin 
thicknesses for asymmetric nuclei. However, in the ERMF model, the mixed interactions of 
a and u mesons with the p mesons makes it possible to give a wide spread in the values of 

nn 

neutron-skin thicknesses without affecting the other properties of finite nuclei |28L 1291] . 

VII. NEUTRON STARS 

In Sec. El we have seen that the energy per nucleon decreases as the OMSC parameter ( 
increases (see FigJJ). This effect is more pronounced at higher densities which implies that 
as ( increases the EOS at higher densities becomes softer. Thus, it is natural to expect the 
differences in the properties of the neutron stars obtained for our new parameter sets with 
different values of OMSC. In this section we present our results for the properties of the 
neutron stars obtained for the seven different parameter sets FSUGZ00 — FSUGZ06. The 
neutron star properties we have considered are the limiting mass, central density and radius 
for the neutron star with the "canonical" mass of 1.4M . These properties are determined 
by integrating the Tolman-Oppenheimer-Volkoff (TOV) equations [30 1 . The TOV equations 
are solved using EOS for the matter consisting of neutrons, protons, electrons and muons. 
The composition of matter at any density is so determined that charge neutrality and (3— 
equilibrium conditions are satisfied. For densities higher than 0.08 fm -3 , the nuclear part of 
EOS is evaluated within the ERMF model and that for electrons and muons the Fermi gas 
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approximation is used. At densities lower than 0.08 fm -3 down to 6.0 x 10~ 12 fm~ 3 we use 
the EOS of Baym-Pethick-Sutherland 31 1. 

In Fig. El we plot mass-radius relationship for the neutron stars calculated using 
FSUGZ00, FSUGZ03 and FSUGZ06 parametrizations. For the sake of comparison we 
also plot the similar results obtained using FSUGold and NL3 parameter sets. The re- 
gion bounded by R ^ 3GM/c 2 is excluded by the causality limit j^. The line labeled by 
AI/I = 0.014 is radius limit estimated by Vela pulsar glitches 33J. The rotation constraint 
as indicated in Fig. |H]is obtained using 34j , 

/ M \ 1/2 /i0km\ 3/2 TT 
, fc = 1045 (-j (_j Hz (25) 



where, the frequency is taken to be 641 Hz which is the highest observed spin rate from the 
pulsar PSR B1937+21. We see from Fig. [HI that as the mass of the neutron stars increases, 
the radius show stronger ( dependence. In particular, radius for the 1.4M Q neutron stars 
decreases linearly from 13.2 Km at £ = 0.0 to 12.4 Km at ( = 0.06. We also calculate the 
radiation radius, 

Roc = ; = (26) 

2GM 



Rc 2 

for the neutron with the canonical mass 1.4M . For our parameter sets, Roo lies in the range 
of 15.19 - 15.93 km. 

In FigElwe have plotted the variation of the limiting mass M max for the neutron stars as 
a function of OMSC parameter (. The ( dependence of the various neutron star properties 
considered here are obtained using the parameter sets FSUGZ00 — FSUGZ06 for which £ 
takes the values 0.0 — 0.06, respectively. The remaining parameters for the FSUGZ sets 
are adjusted to fit the properties of the finite nuclei. The symbols open circle and square 
at ( = 0.0 and 0.06 represent the values of the limiting mass for the NL3 and FSUGold 
parametrizations, respectively. The value of M max decreases with increase in (, because, with 
the increase in ( the EOS becomes softer. It is to be noted that the NL3 and FSUGZ00 
parametrizations have ( = 0. The only difference is that the FSUGZ00 parameter set 
contains additional contributions from the mixed interactions between a, u and p mesons 
and values of K, J and L are much smaller. The value of M max for FSUGold and FSUGZ06 
parameter sets are vary close, since, both these parameter sets have C = 0.06 and have 
similar nuclear matter properties. For our case the M max varies between 2.3M Q — 1.7M Q 
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10] 



for ^ £ ^ 0.06. For the similar range of £, the calculations performed in Ref.[6] yield the 
M max varying between 2.8M — 1.8M Q which is significantly larger. The value of M max in 
Ref. 61 is obtained without including the contributions from the mixed interactions. Unlike, 

n 

in our case, the parameter sets for different values of OMSC used in Ref. |6J are generated 
without constraining the value of L as given by Eq . (|24j) . We must also point out that for 
moderate values of OMSC (( ~ 0.03) we get M max ~ 1.9 M Q . In the present work we do 



not consider the influence of the sel: 
neutron stars. It is shown in Ref. 



-coupling of the p mesons on the maximum mass of the 
6] that when the p meson self-coupling is varied within 
the bounds of the "naturalness", the maximum mass of the neutron star composed of the 
/^-equilibrated matter can increase at most by O.1M . 



In Fig. |H]we plot our results for the proton fractions as a function o: 
solid circles represent the threshold density at which the condition 



baryon density. The 



F n 1/3 ^ F p 1/3 + F e 1/3 (27) 



for the direct Urea (DU) process to occur is satisfied, where, Yi = pi/p is the fraction for 
the i—th species. For the neutron star composed of neutrons, protons, electrons and muons, 
the critical value of the proton fraction Ydu a t which DU process sets in can be obtained 
from Eq. (|2*7jl together with the charge neutrality condition as, 

Y °° = ~ " iTiT- (28) 



1 + 



1+' * 



It is clear from the Eq. (|28|) that at very low densities when = 0, Ydu becomes 0.11. At 
extremely high densities when Y e w = |, Ydu becomes 0.15. From Fig. [Sjwe see that 
the critical proton fraction is ~ 0.13. In FigOwe plot the variations of the densities pu and 
p c as a function of (, where, the pu is the density at which the direct Urea process sets in 
and p c denotes the central density for the neutron star with the "canonical" mass of 1.4M Q . 
As expected, pu and p c increase with (. However, it is interesting to note that pu > p c for 
smaller ( and pu < p c for ( > 0.02. This means that the direct Urea process can occur in 
the neutron stars for the "canonical" mass of 1.4M Q only for ( > 0.02. Also, one can infer 
from the Fig El that the minimum mass Mu for the neutron star in which direct Urea process 
can occur will decrease with In Fig. El we plot Mu versus (. 

We have not considered the hyperonic degrees of freedom in our calculations. The EOS 
with nucleons, hyperons and leptons are softer than the one without hyperons 3(|. The 
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stiffer is the EOS without hyperons, more is the softening effect when hyperons are included 



37 



381 ] . Thus, the results obtained with smaller values of £ will be affected more than 



the ones obtained for larger values of (. As a result, we expect that the variations in the 
maximum mass with ( would certainly decrease. Further, in the presence of hyperons the 
threshold density for the direct Urea process will be pushed up. 
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VIII. CONCLUSIONS 



The effects of OMSC on the properties of finite nuclei and neutron stars are investigated 
using a ERMF model which includes the contributions from all possible mixed interactions 
between the a, u and p mesons upto the quartic order. The simulated annealing method is 
implemented for minimizing the x 2 function required to determine the best fit parameters. 
For a realistic investigation, we generate seven different parameter sets named FSUGZOO, 
FSUGZ01,...., FSUGZ06 for different values of OMSC by adjusting the remaining parameters 
of the model to fit the properties of the finite nuclei. The properties of finite nuclei used 
in the fits are the total binding energies and charge rms radii for several closed shell nuclei 
ranging from normal to exotic ones. In addition, we also include in our fit neutron-skin 
thickness for 208 Pb nucleus to constrain the density dependence of the symmetry energy 
coefficient. The value of neutron-skin thickness is so chosen that it agrees reasonably well 
with the recent predictions All these parameter sets fit equally well the finite 

nuclear properties. But, only moderate values of OMSC are favored from the "naturalness" 
point of view. For our parameter sets the rms errors in the total binding energies calculated 
for the nuclei used in the fits are 1.5 — 1.8 MeV which are significantly smaller than 3.6 MeV 
obtained for recently proposed FSUGold parametrization. 

The behaviour of EOS at higher densities is predominantly determined by the OMSC 
which is different for our different parameter sets. The parameter sets with moderate values 
of the OMSC not only exhibit the "naturalness", but, also yield the EOS for the symmetric 
nuclear and pure neutron matters which closely resemble the ones calculated within the 
framework of the Dirac-Brueckner- Hartree-Fock. For such parameter sets the values of the 
limiting mass for the neutron stars is ~ 1.9M Q . For our various parametrizations the value 
of the limiting mass of the neutrons stars lie in the range of 1.7M Q — 2.3M . We also find 
that the direct Urea process can occur in the neutron stars with "canonical" mass of 1.4M 
only for the moderate and higher values of OMSC. For our various parametrizations, the 
radius of 1.4M neutron stars lie in the range of 12.4 - 13.2 Km. 
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TABLE I: The lower (vo) and upper (vi) limits, maximum displacement (d) and initial values 
(vj n ) for the ERMF model parameters used for minimizing the x 2 value within the SAM. The 
nucleon, omega and rho masses are kept fixed at M = 939 MeV, m u = 782.5 MeV and m p = 770 
MeV, respectively. 





Vn 

u 


Vi 


d 


* in 


e(MeV) 


-16.50 


-15.50 


0.1 


-16.00 


if (MeV) 


220.0 


240.0 


2.0 


230.0 


M*/M 


0.60 


0.70 


0.01 


0.65 


p (fm- 3 ) 


0.145 


0.165 


0.002 


0.155 


9 P 


10.0 


15.0 


0.5 


12.50 


ai(fm _1 ) 


0.0 


0.0095 


0.0010 


0.005 


ai 


0.0 


0.0040 


0.0004 


0.002 


Q2(fm _1 ) 


0.0 


0.0380 


0.0038 


0.019 


a' 2 


0.0 


0.0160 


0.0016 


0.008 


a' 3 


0.0 


0.0160 


0.0016 


0.008 


m CT (MeV) 


490.0 


510.0 


2.0 


500.0 
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TABLE II: Various newly generated parameter sets FSUGZOO, FSUGZ03 and FSUGZ06 for the 
ERMF models. The parameters ~K, a±, a 2 are given in units of fm _1 . The nucleon mass M and 
meson masses m a , and mp are given in units of MeV. The parameters for the NL3 and FSUGold 



sets are taken from Refs. 



Parameters 


NL3 


FSUGold 


FSUGZOO 


FSUGZ03 


FSUGZ06 


tJcr 


10.21743 


10.59217 


10.65616 


10.76145 


11.02412 




12.86764 


14.30207 


13.95799 


14.11104 


14.66595 


yp 


8.94880 


11.76733 


14.32687 


14.67414 


14.52185 


H 


0.019573 


0.007194 


0.030215 


0.015606 


0.006949 


A 


-0.015914 


0.023762 


-0.004544 


0.009753 


0.024487 


c 

s 




0.060000 


0.0 


0.030000 


0.060000 








0.003867 


0.001031 


0.000045 


«'l 


_ 


_ 


0.000779 


0.000507 


0.000053 


a 2 






0.029179 


0.030682 


0.025836 


a' 2 






0.013501 


0.011625 


0.015688 


a' 3 




0.600000 


0.014759 


0.013598 


0.015849 


M 


939 


939 


939 


939 


939 


m a 


508.194 


491.500 


495.763 


500.511 


501.370 


m w 


782.501 


782.500 


782.500 


782.500 


782.500 


m p 


763.0 


763.0 


770.0 


770.0 


770.0 
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TABLE III: The values of parameters expressed as dimensionless ratios using Ea.Hll|). All the 
values have been multiplied by 10 3 





NL3 


FSUGold 


FSUGZOO 


FSUGZ03 


FSUGZ06 


1 

2c 2 , M' 2 
1 

2c' 2 M 2 


1.403 
2.097 


1.221 
1.698 


1.227 
1.782 


1.227 
1.744 


1.173 
1.614 


1 

8c 2 M 2 

K 

6M 


1.030 
0.686 


0.596 
0.252 


0.409 
1.0508 


0.390 
0.547 


0.399 
0.243 


A 
24 


-0.663 


0.990 


-0.189 


0.406 


1.020 


c 

24 




2.500 




1.250 


2.500 


ai 
M 






0.813 


0.217 


0.009 


K 

2 






0.389 


0.253 


0.026 


at2 
4M 






1.533 


1.612 


1.357 


% 
8 






1.688 


1.453 


1.961 


8 




7.500 


1.844 


1.699 


1.981 



TABLE IV: The Nuclear matter properties at the saturation density for newly generated parameter 
sets FSUGZOO, FSUGZ03 and FSUGZ06 are compared with the corresponding ones obtained using 
NL3 and FSUGold parameter sets. The quantities given below are: e the binding energy per 
nucleon, K the nuclear matter incompressibility coefficient, J the symmetry energy, L = 3p^ 
related to the slope of the symmetry energy M* /M is the ratio of the effective nucleon mass to 
the nucleon mass and po the saturation density. 





NL3 


FSUGold 


FSUGZOO 


FSUGZ03 


FSUGZ06 


e(MeV) 


-16.25 


-16.29 


-16.03 


-16.07 


-16.05 


K(MeV) 


271.5 


230.0 


240.0 


232.5 


224.9 


J (MeV) 


37.40 


32.59 


31.43 


31.55 


31.17 


L (MeV) 


118.56 


60.56 


62.19 


64.01 


62.43 


M*/M 


0.595 


0.609 


0.605 


0.603 


0.607 


po(fm~ 3 ) 


0.148 


0.148 


0.149 


0.147 


0.146 
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FIG. 1: The difference between energy per nucleon Ae = e(£ = 0) — e(£) are plotted as a function 
of £ for the SNM (dotted lines) and PNM (solid lines). The results for the e(£) are obtained for 
£ = — 0.06 by using the parameter sets FSUGZ00 — FSUGZ06. The differences Ae at higher 
densities are very sensitive to the value of £. 

FIG. 2: The EOS for the (a) SNM and (b) PNM for the FSUGZ03 parametrization are compared 
with the ones calculated within the DBHF framework for a realistic potential 2^]. Similar results 
for the NL3 and FSUGold parametrizations are also shown. 

FIG. 3: Relative errors in the total binding energy SB = (B ex P - B th )/B ex P for the newly 
generated parameter set FSUGZ03. For the sake of comparison, the values of 5B obtained for 
NL3 and FSUGold parameter sets are also displayed. The rms errors in the total binding energies 
obtained by considering the nuclei used in our fits are 1.6, 2.5 and 3.6 MeV for the FSUGZ03, NL3 
and FSUGold parameter sets, respectively. 

FIG. 4: Relative errors in the charge rms radii 5r c h for the newly generated parameter set 
FSUGZ03. Similar results obtained using NL3 and FSUGold parameter sets are also plotted. 
The rms errors in the charge rms radii obtained by considering the nuclei used in our fits are 0.03, 
0.02 and 0.03 fm for the FSUGZ03, NL3 and FSUGold parameter sets, respectively. 

FIG. 5: Comparison of the results for the neutron-skin thickness Ar = r n — r p obtained using 
the FSUGZ03 , NL3 and FSUGold parameter sets. The crosses are the recent predictions for Ar 



based upon the DBHF calculations 



251 ] and isospin diffusion data 



FIG. 6: Relation between the neutron star mass and its radius R for the FSUGZ00, 
FSUGZ03,FSUGZ06, FSUGold and NL3 parametrizations. The various constraints as indicated 
by causality, rotation and AI/I = 0.014 are discussed in the text. 

FIG. 7: Variation of the limiting mass M max of the neutron star with the OMSC parameter £. 
For different values of £ the remaining parameters of the ERMF model are adjusted to give best 
fit to the properties of the finite nuclei. 
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FIG. 8: Proton fractions as a function of density for the FSUGZOO, FSUGZ03, FSUGZ06, FSUGold 
and NL3 parametrizations. The solid circles indicate the threshold density for the direct Urea 
process. 

FIG. 9: The density py at which direct Urea process sets in and the central density p c for the 
neutron stars with 1AM & as a function of OMSC parameter £. 

FIG. 10: The minimum value of the neutron star mass M\j in which direct Urea process occurs 
as a function OMSC parameter £. 
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Fig. 2b 



1.5 - 



to 



0.5 - 







-0.25 



o 



8 



□ 



o d 
□ o 



o 



6 



o 



8: 



T 



O FSUGZ03 
□ FSUGold 
O NL3 



o 

❖ 



1 



1 







40 



80 



120 160 



8 
o 



200 



Fig. 3 





i 

□ 


1 i 

□ 


i 1 


1 


1 

o 


1 1 1 

FSUGZ03 




1.5 




o 






□o 


FSUGold 
NL3 


- 


1 




□ 










o 


0.5 




e 




□ 

o 






□ ~ 
o 





o 




1 


o 

1 


1 


1 , 1 





40 80 120 160 

A 



Fig.4 



0.6 - 



§ 

o 



o 



o 



□ 

o 



o 



1 



— i 1 r 

O FSUGZ03 

□ FSUGold 

O NL3 

X EXP 



O 



o 



1 



1 



o 



40 80 120 160 200 



Fig. 5 




Fig. 6 



3 




1.5 







0.02 



0.04 



0.06 



Fig. 7 




0.3 0.6 

P [fm 3 ] 




Fig. 10 



